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Abstract 



The N-vortex solutions to the two-dimensional Ginzburg - Landau equa- 
tions for the k = l/v2 parameter are built. The exact solutions are derived 
for the vortices with large numbers of the magnetic flux quanta. The size of 
vortex core is supposed to be much greater than the magnetic field penetration 
depth. In this limiting case the problem is reduced to the determination of 
vortex core shape. The corresponding nonlinear boundary problem is solved 

by means of the methods of the theory of analytic functions. 
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I. INTRODUCTION 



In different areas of physics there is a wide interest in nonlinear field equations admitting 
particlelike solutions. The Ginzburg-Landau (GL) theory is a well known example of this 
kind. In the case of a hard superconductor subjected to an external magnetic field the GL 
equations give vortex-line solutions. The properties of the vortex states in the GL theory are 
investigated in two limiting cases |2| . In the first case the sector of the fields near the lower 
critical field H Cl is considered. In this sector a mixed state of the type-II superconductor is 
a system of widely separated vortices. In this case the problem is solved for the unit vortex. 
Due to the cylindrical symmetry of the system the corresponding GL equations are reduced 
to the set of the ordinary differential equations. In the second case the sector of the fields 
near the upper critical field H C2 is considered. In this region the periodical lattice of the 
vortices is the equilibrium state of the superconductor. Due to the small value of the order 
parameter the GL equations turn out to be linear. It allows one to investigate them by the 
usual methods. 

In this work the exact solutions of nonlinear GL equations are built. The solutions are 
derived in the special case of the k = S/£ = l/y/2, where S is the magnetic field penetration 
depth, £ is the fluctuation radius of the order parameter, and in the asymptotic of the large 
magnetic fluxes of each vortex $ 3> $o? where $o — Tthc/e is the magnetic flux quantum. 
In the limiting case when the vortex core sizes considerably exceed the typical lengths 5 or 
£ the problem is reduced to the purely geometrical one, to the determination of the shape of 
the contour confining the vortex cores. Such a problem proves to be nonlinear. Nevertheless, 
it is solved by means of methods of the theory of analytic functions. 

This case is interesting since the quantum properties of the superconductor are mani- 
fested in the classical limit. The k parameter determines the value of tension coefficient on 
the interphase of normal and superconducting phase. It is at k = 1/ a/2 that the coefficient 
vanishes |J and causes the agreement of all three critical magnetic fields H Cl = H C2 = H c . 
Therefore, the material with such value of the k parameter possesses the properties of both 



2 



type-I and type-II superconductors. On the one hand, the vortices with large magnetic fluxes 
are the incorporations of the normal phase. They should be considered as the parts of the 
intermediate state of type-I superconductor. On the other hand, such regions of the normal 
phase are nothing but the vortices. As shown below their shape depends on the coordinates 
of the zeros of the order parameter and the numbers of the magnetic flux quanta in every 
vortex, i.e. on the typical parameters of the mixed state. 

The vortex states with large magnetic fluxes under consideration are not too "exotic" 
objects. Parallel with the macroscopic vortices the quantum vortices can also exist in su- 
perconductor. This is caused by the absence of the interaction among the vortices. The 
energy of the whole system depends only on the total magnetic flux ||. At the given total 
magnetic flux the decay of the large vortices into small vortices does not lead to a gain in 
energy. 

This consideration is based on the results of Ref . . In this work it is shown that the GL 
equations are reduced to a set of the first order coupled equations for the vector potential 
and the order parameter. By appropriate choice of the variables the latter are written in 
the form of an equation similar to the Poisson equation [[5] . Its form is exactly the one that 
has the asymptotically exact solutions ||. 



Let us consider the two-dimensional distributions of the complex order parameter ip(r) 
and the vector potential A(r). Such distributions are realized for the superconducting 
cylinder subjected to the longitudinal magnetic field (see Fig. |l|). If no edge effects are 
present all, values depend on the transversal coordinates x and y only. 

After rescaling the lengths and the fields the free energy of superconductor per unit 
length of the sample is written as follows 



II. GENERAL APPROACH 
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(1) 



where all the lengths are measured in the units of magnetic penetration depth 5, and A = \^2k 
is the dimensionless coupling constant. For the chosen variables the physical magnetic field 
B = H c corresponds to the dimensionless field B = A/2. 

It is convenient to introduce the complex variables such as 

z = x + iy, z = x — iy (2) 

and the vector potential 

A=- (A x - iAy) , A = - (A x + iAy) . (3) 

These notations us permit to present the equations derived in Ref. [|,|J for k = l/y/2 and, 
correspondingly, A = 1 in the most simple form. Such equations give the minimum of free 
energy ([]]). Unlike the GL equations they are the first order with respect the derivatives: 

H - iA^ = 0, (4) 
1 /, M 2 



s + -ij = o. (5) 

These equations give the vortexlike solutions with the magnetic flux 

$ = J d 2 rB = 2vrn, (6) 



where n is the positive integer number. The corresponding total free energy of the system 
in dimensionless units turns out to be equal exactly to irn. 

The set of the coupled equations (|j) and @ is reduced to one equation by the substitution 
of the form : 

A = -l^(« + 2]n|/|), (7) 
^ = ^exp(-^), (8) 
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where u(x, y) is the real- valued function, and f(z) is the analytic function. After the substi- 
tution the former equation of the system becomes the identity. The second equation takes 
the form 

Au = -4vrp (u) — 2A In |/| 

N 

= -4irp (u) - 47r^n j (5(r - r^), (9) 

3=1 

where A is the Laplace two-dimensional operator, 

p{u) = i- [exp(-u) - 1] . (10) 

Formula @ differs from the analogous one of Ref. || in the last term of the right-hand 
side. Here singularities of the function In \f(z) | are taken into account. The analytic function 
f(z) vanishes like (z — Zj) Uj at the positions of the vortices Zj, where nj is the positive integer 
number, j = 1, N, and N is the number of the vortices 0,0]. This leads to the S functions 
in the right-hand side of the equation ([|). 

The boundary condition for this equation is vanishing of the value u(x,y) at the infinity. 
According to expression (||) this requirement gives equality = 1. The function f(z) is 
to some degree an arbitrary one. Its form depends on the choice of the gauge of the vector 
potential. 

The derived equations have the following physical meaning. They evidently describe the 
screening of the external positive point charges in a classical plasma. The value of the u(r) 
is the electrostatic potential of the system of charges. The dependence of p(u) gives the 
induced charge density described by the jellium model which consists of a negative smeared 
background and mobil charged fluid. 

The exact solutions of nonlinear screening equation (|9|) are built in the limiting case 
rij ^> 1. In this approximation a sufficiently large positive potential is supposed to exist near 
external point charges. As a consequence the screening results from a complete repulsion 
of the positively charged fluid from the regions surrounding the point charges. The induced 
charge density turns out to be approximately constant within these regions and is equal to 
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p(oo) = —1/Att. The size of such region is much greater than the Debye length for linear 
screening. 

In these assumptions one can ignore the details of the potential behavior on the scale of 
the order of the Debye length. Namely, one can suppose that the system consists only of 
the uniformly charged spatial region and a set of point charges. In this case the screening of 
the point charges can be achieved by means of selecting such a form of the region in which 
everywhere outside the region both the potential and the electric field would vanish. This 
problem is much easier than the initial one. Due to its two- dimensionality it is solved by 
the methods of the theory of functions of a complex variable ||. 

In the given case the results of the work || can be rewritten in the following form. The 
solution is given by the function 

^(0 = co + E (ii) 

j=l S S>j 

which makes a conformal mapping of the inner region of the unit circle |£| < 1 in an 
auxiliary plane of the complex variable ( into the region of interest. The constants Q obey 
the condition \Q\ > 1 for all j. The shape of the contour is specified parametrically: z = u>(() 
for ( = exp(za) and < a < 2n. 

The set of unknown constants Cj and (j is calculated by means of the system of nonlinear 
algebraic equations 

u{Q=z v (12) 

where £ 3 - = l/Q-, j = l>--->N; oj' (C) is the derivative of the function uj((). Constant Co is 
determined from the additional requirement that an inner point of the unit circle should 
correspond to a given point within the uniformly charged region. 
The potential u (x, y) is given by the expression 

u(x,y) = ^\z\ 2 -2ReW(z) } (13) 
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where 

1 N 

W(z) = -V(z) + ^2 n 3 m 
8 3=1 

Here function C{z) is the inverse of z — 

The results of the electrostatic problem allow one to write the corresponding solutions of 
the GL equations easily. In the next section the general formulas are given and some simple 
cases are considered. 



QC (g) - 1 

C (*) - 



(14) 



III. RESULTS AND DISCUSSION 



Taking into account the fact that the gauge of the vector potential can be an arbitrary 
one, let us take the function f(z) in the simpliest form 
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m = n 



Z — Zn 



(15) 



Within the core of the vortices the order parameter and the vector potential will be described 
by the expressions: 

v r z-Zj U(z)-1 



\Z — Za 



x exp 



\z\ 2 — Ret> (z) 



A= — 



N 



n.; 



j=l Z Z 3 



(16) 
(17) 



where 
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In formula (P^7| ) value A will be regular near the points z — Zj, because the singularities of 
the function r](z) are eliminated by the third term in square brackets. 

Outside the vortex cores the order parameter and the vector potential are described by 
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a = — 2 Y,— '-■ (20) 

In this region the induction of magnetic field B equals zero. On the contrary, within the 
vortex core the induction is B = 1/2. The last value corresponds to the field B = H c 
measured in the physical units. 

Formulas (|T6|)-(pO|) give the solutions only when the vortex cores are connected to each 
other U . For the sets of widely separated vortices these formulas should be applied to each 
connected group. These expressions do not work for the regions of "Swiss cheese" type. A 
solution of an electrostatic problem of this kind was obtained in Ref. ||. 

Formulas (|T6|) -(pO|) present the desired results as the functions of the auxiliary variable 
£. To obtain the explicit dependence of these functions on the variable z it is necessary 
to invert the equality z = ou(C) with respect to the variable (. It is possible only for some 
vortices. 

Let us consider two special cases. The simpliest one is the problem for the unit vortex 
in the origin of the coordinates z = 0. The corresponding solutions have the form 







(;)" exp 









(21) 



A = —-z. 



where \z\ < r. The vortex core is a circle with the radius r = 2 v / n, where n is the number 
of the flux quanta of the magnetic field. 

The problem for the twin vortices placed at the points x = ±a of the real axes can be 
regarded as the second example. The vortex cores are overlapped if the distance a obeys 
the inequality a < r. Their shape is described by the conformal mapping function 

where 
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ci = a^(p 2 - 1) [p+^p 2 - lj. 



Ci 
p 



p+\/p 2 - 1, 



(23) 



The order parameter is given by the expression 



+ V* 



a 2 + r z 



x exp 



+ -Re ( V 



(24) 



The corresponding formula for the vector potential is too bulky to be presented here. 

Note that solutions ([|J)-(|l7D may be rewritten in a form similar to the eigenfunction of 
a charged particle in the lowest Landau level: 

/ i 1 2 \ 



ip = F(z) exp 



(25) 



where F(z) is an analytic function. Such a presentation is evident for the following physical 
reasons. Within the normal phase the magnetic-field induction is approximately constant 
and equal to the critical field H c . The order parameter is negligibly small here as well. The 
GL equations linearized with respect to the smallness l^l <C 1 really give such solutions. 
To derive the presentation ( p5|) let us take function f(z) in the form f(z) = exp (W(z)). 



Correspondingly, formula (|i~3"D gives the equality F(z) = f(z). The vector potential takes 
the form A = —iz/8. 

In conclusion, let us discuss the validity of the model. It was assumed above that the 
order parameter depends on two coordinates only. This distribution is not realized in fact. 
According to the theory of the intermediate state of superconductors the domains of a 
normal phase decay near the surface of a sample into the alternating layers of normal and 
superconducting phases 0. Thus, there is a dependence on the third coordinate directed 
along the external magnetic field. In this connection the existence of the macroscopic vortices 
may be provided only by a special mechanism for the attraction of vortices. The pinning 



of vortices near space inhomogeneities or free boundary of the sample can become such a 
mechanism. The solution of this problem is beyond the scope of this article. 
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FIGURES 

FIG. 1. The geometry of the system. The shaded region of the picture shows a unit vortex in 
a cylindrical specimen of superconductor. 
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